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Abstract 

A new scheme to calculate the exchange tensor J_^. describing in a phenomenological way the 
anisotropic exchange coupling of two moments in a magnetically ordered system is presented. 
The ab-initio approach is based on spin-polarised relativistic multiple-scattering theory within 
the framework of spin-density functional theory. The scheme is applied to ferromagnetic CrTe as 
well as the diluted magnetic semiconductor (DMS) system Gai_^Mn^As. In the later case the 
results show that there is a noticeable anisotropy in the exchange coupling present, although not 
as pronounced as suggested in recent theoretical investigations. 

PACS numbers: Valid PACS appear here 
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I. INTRODUCTION 



The mapping of the energy of a magnetic sohd calculated from first principles for different 
magnetic configurations onto a Heisenberg Hamiltonian is nowadays a widely used concept 
that allows a number of interesting subsequent investigations. An example for this is the de- 
termination of the Curie temperature of a ferromagnet by means of Monte Carlo simulations 
using the calculated exchange coupling parameters J^- as input (z and j are indices labeling 
the individual atomic sites) [l]. Apart from fitting the Ji/s to the total energies obtained 
for different magnetic configurations one can use the energies of spin spirals as basis for such 
a mapping As an alternative one may use perturbation theory, that allow to calculate 
the Jjj's directly. In fact the expression derived by Lichtenstein et al. [sl within the frame 
work of non-relativistic multiple scattering theory is now successfully used for a wide range 
of materials 

Initiated among others by investigations on the magnetic ground state configuration of 
nano-scale systems, there is strongly growing interest in the interplay of exchange coupling 
and spin-orbit coupling [7l|. Besides the magnetic anisotropy energy, the spin-orbit coupling 
gives rise to an anisotropic exchange coupling. Using again the above mentioned concept 
the isotropic exchange constants Jij have to be replaced by a corresponding exchange cou- 
pling tensor J. .. By generalising the approach of Lichtenstein et al. to a fully relativistic 
formulation, Udvardi et al. [8] could derive corresponding expressions for the elements of 



J... Corresponding applications to thin films as well as to finite deposited clusters can be 
found in the literature [8, I^. A disadvantage of the expressions worked out by Udvardi et al. 

is that one has to use various magnetic configurations as a reference state to determine 
all elements of J_._.. In the following we present an alternative approach that can be derived 
in a rather transparent way and does not have that problem. As it is demonstrated both 
approaches give nevertheless results that are quite close to each other. As an application 
of our new scheme we present results for the exchange tensor J_._. in ferromagnetic CrTe 
and in the diluted magnetic semiconductor (DMS) system Gai _xMn,: As. The later will be 



discussed in relation to recent work of Timm and MacDonald [10|] who used a tight-binding 
description of the system in contrast to the ab-initio approach employed here that is based 
on local spin density approximation (LSDA). 
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II. THEORETICAL APPROACH 



Starting point of our derivation for J__ is the expression for the change in energy lS.Eij 
of a system upon a perturbation taking place at sites i and j. An expression for AiJ^ was 
worked out by several authors [sl, [ll| within the framework of multiple scattering theory 
and making use of Lloyd's formula. The derivation of this expression can straight forwardly 
be applied when working in the framework of spin-polarised relativistic multiple scattering 

n 

or KKR [12] formalism. Adopting the convention for the corresponding electronic Green's 



function as used by Dederichs and coworkers 



131 ] its off-site part is given by: 



G(rl,r,,E) = -ipY,R\{n,E)Gl^,{E)R^^{r,,E) , (1) 

AA' 

where G^l^,{E) is the so-called structural Green's function, R\ is a regular solution to the 
single-site Dirac equation labeled by the combined quantum numbers A (A = with 
K and being the spin-orbit and magnetic quantum numbers |l4l) and p is the electron 
momentum. The energy change AEij can then be written as [3], lll| 

AEij = -i$5 j dETiaceAfG'^APG^' , (2) 

where Af is the change of the single-site t-matrix due to the perturbation Ay*(r) at site i 
and the underline denotes matrices with respect to the quantum numbers A. To first order 
in AV^{r) the change Af is given by 

At\,^ = J d'rR^,{r)AV{r)Rl{r) = AVifl' . (3) 

Using instead the convention for the Green's function as used by Gyorffy and coworkers 
[isl ] one may express AEij in terms of the scattering path operator t^I,j^{E) 

AEij = j dETi&ceAV^^^Y^ AV}^^^T^' , (4) 

where use have been made of the relation G^^ = {t'')~^T^^ (P)'^ for i ^ j and the matrix 
(^)* "--'"luated using the alternative set of regular solutions Z\ to the 



elements AV^^; are to be eva. 



single-site Dirac equation 
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Changing the orientation of the spin magnetic moment fhi within an atomic cell i and 
adopting the rigid spin approximation (RSA) 1^ implies a corresponding change of the 
spin-dependent potential (3aB{r) by: 



^V{r) = Vf.{r) - Vf,,{r) = (3a{n - ho)B{r) , (5) 

where P is one of the standard Dirac matrices and a is the vector of 4 x 4-spin matrices 
[igI . In writing Eq. ([5]) a collinear spin magnetisation within the cell has been assumed 
together with a change of its orientation from ho to h. Accordingly, B{r) corresponds to 
the difference of the spin-projected potential functions B{r) = ^{V^{r) — V^{r)) 12|. This 
leads for the matrix elements AV^^)* to: 



^vS' = E ^yiT^o. (6) 

with 



a=x,y,z 



Ayjf)- = I dhZl{r)(3a^B{r)Zr,,{r) . (7) 
Comparison with the generalised Heisenberg Hamiltonian 



Hex = (8) 

with ej(j) the orientation of the spin magnetic moment at site allows one to write for 
the elements of the exchange coupling tensor J_.. 

J";"' = --53 [ dETraceAV^^^''^T'^AV^^^''^r'' . (9) 

^ TT J 

The scheme outlined above has been implemented using the spin-polarised relativistic 
(SPR) version of multiple scattering theory AH calculations have been done within 

the framework of the local spin density approximation (LSDA) to spin density functional 



theory 17|]. To represent the results for the exchange tensor J__ we use_ the conventional 
decomposition of the corresponding Heisenberg Hamiltonian in Eq. ([8]) 
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FIG. 1: Isotropic exchange interaction parameters Jij between Cr atoms at sites i and j in ferro- 
magnetic CrTe as a function of the inter-atomic distance Rij. The results based on the present 
approach (full symbols) are compared to results obtained using the approach of Udvardi et al. ^ 
(open symbols). Circles and squares represent the coupling of a Cr atom in layer 1 (Crl) to another 
Cr atom in layer 1 (Crl) or layer 2 (Cr2), respectively. 



ij ij 

Here Jij is the isotropic exchange coupling constant, J_^. is the traceless symmetric part of 

— '^j 

J_. . and the antisymmetric part is represented by the Dzyaloshinski-Moriya (DM) vector Dij. 
— u 

It should be emphasized that isotropic in the context of Jij refers to spin-space and does not 
imply that there in no anisotropy in real space, i. e. Jij will in general not only depend on 
the distance \Rij\ between two sites but also on the orientation Rij of the distance vector. 

III. RESULTS AND DISCUSSION 

To demonstrate the application of our approach we present in Figs. [U and [2] results for 
the coupling parameters Jij and Dij of the two inequivalent Cr atoms in ferromagnetic CrTe 
having the NiAs structure. The isotropic parameters Jij shown in Fig. [1] reflect dominating 
ferromagnetic coupling that is quite far reaching, i.e. slowly decaying. As one notes the 
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FIG. 2: Components of the Dzyaloshinski-Moriya interaction vector Dij (from top to bottom: 
X, y, z) between Cr atoms at sites i and j in ferromagnetic CrTe as a function of the inter-atomic 
distance Rij. The results based on the present approach (full circles) are compared to results 
obtained using the approach of Udvardi et al. (open squares). 
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FIG. 3: The same as in Fig. [2] top but with spin-orbit couphng scaled by factor 2. 



isotropic exchange couphng between a central Cr atom in layer 1 (denoted Crl) to another 

Cr atom in layer 1 and 2 (denoted Crl and Cr2 and represented by squares and circles, 

respectively, in Fig. [T]) is in the same order of magnitude. This means there is no remarkable 

spatial anisotropy imposed by the layered structure of the system for the isotropic coupling 

constant Jij. The anisotropy of the exchange coupling is represented by J_^. as well as by Dij. 

— 

As J_^. turns out to be quite small we show in Fig. [2] only the three components of the DM 
— 

vector Dij. Many of the DM vector components are zero due to symmetry restrictions 18|. 
In particular one finds a non-zero DM vector Dij only if the sites i and j belong to different 
sub- lattices Crl and Cr2. Due to the NiAs structure of the system the non- vanishing x- 
and y-components of the vector are of the same order of magnitude while the z-component 
is one order of magnitude larger. The different behaviour of the x-, y- and 2;-components 
reflects obviously to some extent the quasi-layered structure of the system (hexagonal Cr 
layers with Te-layers in between). However, the anisotropic exchange coupling is still about 
two orders of magnitude smaller than the isotropic one. 

As the comparison of the results for Jij and Dij obtained using the approach presented 
above and that of Udvardi et al. [8|], respectively, in Figs. [T] and [2] demonstrates, both 
schemes give very similar results. This also holds for other systems studied so far with most 
pronounced differences occurring for the DM vector. It should be stressed, however, that the 
above scheme allows to determine J with respect to one common reference state; i.e. there 
is no need to use various reference states to get all tensor elements. This ensures that the 
results for the various elements are always consistent even when the choice of the reference 
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FIG. 4: Isotropic exchange interaction Jij between Mn atoms at sites i and j in Gai-^^Mn^-As, 
scaled by the factor (Rij/a)'^, as a function of the inter-atomic distance Rij. Results are given for 
7 at.% Mn (full circles) and for 4 at.% Mn (open squares). 



state influence the result, e.g. when the RSA is not fully justified. 

To demonstrate that the DM interaction is indeed induced by spin-orbit coupling (SOC) 
we performed model calculations with the strength of SOC artificially increased by a factor 
of 2. While the isotropic exchange coupling constants Jij hardly changed, the DM vector 
components increased essentially by the same factor. This can be seen by comparing Df^ 
given in Fig. [3] with the results in Fig. [2] (top). Fig. [3] also shows that for this specific 
situation results based on the two approaches considered may differ in an appreciable way. 

Figs, m and [5] show the results of an application of our approach to the diluted magnetic 
semiconductor system Gai_a;Mn^As. The isotropic exchange coupling (Fig. H]) and also its 
concentration dependence agrees quite well with the results of other authors jq] , indicating in 
particular that the spin-orbit coupling accounted within the present work affects the isotropic 
exchange coupling only slightly. As mentioned above and as was noted by other authors 
[gI. [lol there is a directional dependency for Jij. This is demonstrated by Fig. [6]where results 
for the concentration x = 0.04 are given for R^j along [001], [110], and [111] separately. As 
for the magnitude of Jij, this spatial anisotropy of J^j is only slightly influenced by inclusion 
of SOC. In contrast to this SOC is ultimately responsible for the anisotropy in the exchange 
coupling represented by the DM vector shown in Fig. [5l In contrast to the CrTe system 
considered above the three components of the DM vectors are of the same order of magnitude 
as a consequence of the zincblende structure of Gai_a;Mna,As. 
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Although, it seems not possible to give a simple scaling behaviour of the magnitude of 
the exchange coupling parameters with respect to the inter-atomic distance Rij, one notes 
that the components of Dij decay less rapidly as Jij with increasing Rij. This behaviour was 
also found for other systems and is in line with the findings of Timm and MacDonald [10|] . 
However, our results for the isotropic as well as anisotropic exchange coupling constants 
differ quite appreciably from those obtained recently in a more phenomenological way by 
these authors. In particular the tensor elements representing anisotropic exchange are found 
to be around one order of magnitude smaller than given in the previous work. The anisotropy 
of the DM interaction is demonstrated in Fig. [7] for the component Dfj. As one notes Dfj 
depends quite strongly in the direction Rij. In particular one finds D^j to be zero e.g. for 
the [111] direction due to symmetry. Also because of the symmetry of the system one finds 
for each direction a symmetry-related one for which the sign of Dfj is reversed. 

The presence of a non-collinear ferromagnetic structure in Gai-a^Mn^^As is assumed to be 
partially responsible for the missing of remanent magnetisation observed experimentally in 
annealed samples {20, 21, 22\- In particular, the remanent magnetisation in this DMS system 
can be noticably increased in the presence of a rather small magnetic field. Such a behaviour 



could indeed be explained by the presence of non-collinear magnetism in the system 



22|. 



The anisotropy o 



authors 



Id, 



23 



the exchange coupling in Gai-ajMn^^As was studied theoretically by various 



24, 



251 ] to find whether it can be responsible for the formation of a non- 



collinear ferromagnetic structure as a ground state in this DMS system. However, these 
investigations were based on phenomenological or semi-phenomenological approaches and 
the results obtained are rather controversial. In contrast to this, the present approach allows 
us to evaluate the elements of the exchange coupling tensor (in particular, its antisymmetric 
part representing the DM coupling) on the basis of ab-initio electronic structure calculations. 
As was demonstrated above, this leads indeed to a rather large value for the DM coupling 
in Gai_a;Mn^.As, which is only about one order of magnitude smaller than for the isotropic 
exchange. As mentioned, this finding is in line with the results of Timm and MacDonald [10]. 
Obviously, the values for the DM coupling term cannot be considered as negligibly small, 
and as a consequence one cannot exclude a noticable non-collinear ferromagnetic order in 
the system. To clarify this question corresponding Monte Carlo simulations based on the 
calculated exchange tensor will be performed. 
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IV. SUMMARY 



A new scheme to calculate the exchange coupling tensor J_^. has been presented that 
is based on ab-initio electronic structure calculations using spin-polarized fully relativistic 
multiple scattering theory and spin density functional theory. Application to ferromagnetic 
CrTe as well as to other systems demonstrates that the approach gives results for the ex- 
change tensor elements very similar to those obtained using the approach of Udvardi et al. 
However, the new approach makes use of a unique reference state ensuring the internal con- 
sistency of the tensor elements. Application to the diluted magnetic semiconductor system 
Gai_a;Mn3;As led to an isotropic exchange in full accordance with previous non-relativistic 
calculations that were also based on ab-initio electronic structure calculations. The results 
obtained for the anisotropic exchange coupling are in accordance with the data of Timm 
and MacDonald concerning the variation with distance. However, numerically the coupling 
constants obtained by the present ab-initio approach and the more phenomenological scheme 
of these authors differ in an appreciable way. 
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FIG. 5: Components of Dzyaloshinski-Moriya interaction vector Dij (from top to bottom: x, y, z) 
between Mn atoms at sites i and j in Gai_j^Mna;As, scaled by the factor (Rij/a)'^, as a function 
of the inter-atomic distance Rij. Results are given for 7 at.% Mn (full circles) and for 4 at.% Mn 
(open squares). 
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FIG. 6: Scaled exchange interactions {Rij/a)'^Jij in ferromagetic Gai-ajMn^^As, x = 0.04, as a 
function of the inter- atomic distance Rij along different directions. 
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FIG. 7: Components of Dzyaloshinski-Moriya interaction Dij between Mn atoms at sites i and j 
in Gai_x.Mn2:As, x = 0.04, scaled by factor (Rij/a)'^, as a function of the inter-atomic distance Rij 
along different directions. 
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